ASYMPTOTIC BEHAVIOUR OF SELF-CONTRACTED PLANAR CURVES 
AND GRADIENT ORBITS OF CONVEX FUNCTIONS 



ARIS DANIILIDIS*, OLIVIER LEY, AND STEPHANE SABOURAU 

Abstract. We hereby introduce and study the notion of self-contracted curves, which encom- 
passes orbits of gradient systems of convex and quasiconvex functions. Our main result shows 
that bounded self-contracted planar curves have a finite length. We also give an example of a 
convex function defined in the plane whose gradient orbits spiral infinitely many times around 
the unique minimum of the function. 
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1. Introduction 

This work is mainly devoted to the study of the length of bounded trajectories of the gradient 
flow of convex (or quasiconvex) functions in the plane. The motivation for this study comes 
from a well-known result due to S. Lojasiewicz (see [E]), asserting that if / : — > M is a 
real-analytic function and x G /~^(0) is a critical point of /, then there exist two constants 
p G [1/2, 1) and C > such that 

||V/(x)||>C|/(x)|'' (1.1) 

for all X belonging to a neighborhood U of x. An immediate by-product is the finite length of 
the orbits of the gradient flow of / lying in U . The proof is straightforward using (jl.ip and is 
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illustrated below: let 7 : [0, +00) — > C/ be a gradient trajectory of /, that is, 7(t) = — V/(7(t)). 
Then, 



The aforementioned inequality (jl.ip has been extended by K. Kurdyka in [13] for functions 
belonging to an arbitrary o-minimal structure (we refer to [11] for the relevant definition), in a 
way that allows us again to deduce the finiteness of the lengths of the gradient orbits in this more 
general context. In [3j and [4j, a further extension has been realized to encompass (nonsmooth) 
functions and orbits of the corresponding subgradient systems. 

It should be noted that in the above cases the functions enjoy an important structural property 
(o-minimality) and that, for general functions, there is no hope to prove a result like p.2|) . A 
classical example of J. Palis and W. De Melo ( |16l page 14]) asserts that the bounded trajectories 
of the gradient flow of an arbitrary C°° function need not converge (in particular, they are of 
infinite length). In the aforementioned example the critical set of the function is not reduced to 
a singleton: in Section 17. H we provide another example of a smooth function having a unique 
critical point towards which all corresponding orbits converge, but again are of infinite length. 

The case when / is a convex coercive function is particularly interesting in view of its potential 
impact in numerical optimization (see [1], [3], [5], for example). But convex functions are far 
from being analytic and they do not satisfy neither the Lojasiewicz inequality nor its generalized 
form established by Kurdyka, unless a growth condition is assumed (see [3 Sections 4.2-4.3] 
for a sufficient condition and a counter-example). Nevertheless, their rigid structure makes it 
natural to think that the orbits of their gradient flow are of finite length. It is rather surprising 
that the answer of this question is not yet known in the literature except in some particular 
cases. 

Let us mention that in the framework of Hilbert spaces, this has been stated as an open 
problem by H. Brezis [6", Open problems, p. 167]. In infinite dimension, R. Bruck [7J proved that 
the (sub)gradient orbits of convex coercive functions are converging towards a global minimizer 
of / but this convergence holds only with respect to the weak topology. Indeed, B. Baillon [2] 
constructed a counterexample of a lower semicontinuous convex function / in a Hilbert space 
whose gradient orbits do not converge for the norm topology. A straightforward consequence 
is that these orbits have infinite length. Concurrently, there are some cases where a convex 
coercive function f : H is known to have (sub)gradient orbits of finite length. This is true 
when the set of minimizers of / has nonempty interior in the Hilbert space H (see H. Brezis [6]), 
or whenever / satisfies a growth condition. For a detailed discussion and the proofs of these 
facts, we refer to [U Section 3-4]. 

The aforementioned results do not cover the simplest case of a convex smooth function defined 
in the plane and having a unique minimum. One of the main results of this work is to prove the 
following: 

Theorem 1.1 (Convex gradient system). Let f : be a smooth convex function with a 

unique minimum. Then, the trajectories 7 of the gradient system 




yielding 





at) 



V/(7(t)) 



have a finite length. 
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The proof of this result does not use the whole convexity of / but, instead, rather relies on the 
convexity of its level-sets. More precisely, the conclusion of Theorem 11.11 remains also true for 
the orbits of the gradient flow of a quasiconvex function (see Corollary I6.3p . 

Actually, both results will follow as consequences of a much more general result (Theorem ll.3p 
about bounded self- contracted planar curves, which allows us to provide a unified framework for 
this study. 

Definition 1.2 (Self-contracted curve). A curve 7 : / ^ defined on an interval I of [0, +00) 
is called self- contracted, if for every ti < ^2 ^ ^3) with tj G /, we have 

dist(7(ti), 7(t3)) > dist(7(t2), 7(i3))- (1-3) 
In other words, for every [a,b] C /, the map t £ [a,b] ^ dist(7(t), 7(6)) is nonincreasing. 

We prove the following. 

Theorem 1.3 (Main result). Every bounded continuous self- contracted planar curve 7 is of 
finite length. More precisely, 

length(7) < (87r + 2)D(7) 
where D{j) is the distance between the endpoints ofj. 

Let us finally mention that, even if gradient orbits of convex functions have finite length in 
the plane, they do not enjoy all the properties of the gradient orbits of real-analytic functions. 
Indeed, on the one hand, the so-called Thom conjecture for the gradient orbits of real-analytic 
functions holds true: ifxoo denotes the limit of the orbit 7(i), then the secants i'jit)—Xoo)/\hit) — 
Xoo\\ converge towards a fixed direction of the unit sphere (see K. Kurdyka, T. Mostowski and 
A. Parusinski [H]). On the other hand, as we show in Section 17.21 an analogous result fails 
in the convex case. Indeed, the orbits of a convex gradient flow may turn around their limit 
infinitely many times. 

Our techniques only work in the two-dimensional case. We do not know whether Theorem ll.il 
and Theorem 11.31 hold in greater dimension. 

The article is organized as follows. In Section [21 we present basic properties of self-contracted 
curves. In Section [3l we decompose each polygonal approximation of a bounded self-contracted 
curve in an annulus centered at its endpoint into horizontal and vertical segments. We establish 
upper bounds on the total length of the vertical segments in Section [3] and on the total length 
of the horizontal segments in Section [H The proof of the main result is presented in Section [5l 
In Section [6l we show that the orbits of various dynamical systems are self-contracted curves. 
Two (counter)-examples are presented in Section [71 

Notations. Throughout the manuscript, we shall deal with the flnite-dimensional Euclidean 
space M" equipped with the canonical scalar product (■,•). We denote by || ■ || (respectively, 
dist(-, •)) the corresponding norm (respectively, distance). Therefore, the distance between two 
points X and y of M? will be denoted by \\x — y\\, dist(x,y) or sometimes \xy\. We also denote 
by dist (x, S) the distance of a given point x G M" to a set 5 C W^, by B{x,r) the closed ball 
with center x G M"' and radius r > and by S{x, r) its boundary, that is, the sphere of the same 
center and the same radius. For < r < i?, we denote by 

U{r, ii) := {x G M" I r < \\x\\ < R} (1.4) 

the annulus centered at the origin O with outer radius R and inner radius r and by AR = R — r 
its width. Let 

\p,q] := {p + t{q-p) |tG [0,1]} 



4 



A. DANIILIDIS, O. LEY, AND S. SABOURAU 



be the closed segment with endpoints g S M". A subset S of M" is cahed convex, if \p,q] C S 
for every p,q & S. 

2. Self-contracted curves 

Throughout this paper, we shall deal with curves 'j : I ^ M", defined on an interval / of M. 
We recall that the length of a curve 7 : / ^ is defined as 

length(7) = sup | dist(7(ti), 7(ti+i))| 

where the supremum is taken over all the finite subdivisions {ti}^^^ of /. 
We shall need the following definition. 

Definition 2.1 (Convergence of a curve). A curve 7 : / — > M*^ is said to converge to a 
point xq € M" if 7(t) converges to xq when t goes to t+ := sup I. 

A curve 7 : / ^ M*^ is said to be bounded, if its image 7(1) is a bounded subset of M". 

We start with an elementary property of self-contracted curves. 

Proposition 2.2 (Existence of left/right limits). Let 7 : / 1-^ M" be a bounded self- contracted 
curve and (a, b) C /. Then, 7 has a limit in when t G (a, b) tends to an endpoint of (a, b). 
In particular, every self- contracted curve can be extended by continuity to the endpoints of I 
{possibly equal to ±00). 

Proof. Since 7 lies in a compact set of M", there exists an increasing sequence {ti} in (a, b) with 
ti ^ b such that j{ti) converges to some point of M", noted 7(6)"'". Fix any i,j € N* and let 
ti < t < ti^j. By (|1.3p . we have 

dist(7(t),7(ti+j)) < dist(7(ti),7(ti+j)). 
Letting j go to infinity, we derive 

dist(7(t),7(fc)+)<dist(7(t.),7(fc)+) 

which gives ^{t) 'jib)'^- 

Further, using the triangle inequality and the inequality (|1.3p . we have 

dist(7(ti), 7(^2)) < dist(7(ti), 7(t3)) + dist(7(t3), 7(^2)) < 2 dist(7(ti), 7(t3)). 
Using this inequality, we can show, as previously, that 7(t) converges as (a, b) 3 t ^ a. 

The last part of the assertion is straightforward. □ 

The following result is a straightforward consequence of Proposition 12. 2i 

Corollary 2.3 (Convergence of bounded self-contracted curves). Every bounded self- contracted 
curve 7 : (0, +00) M" converges to some point xq € as t — > +00. Moreover, the function 
t I— s- dist(xo, 7(t)) is nonincreasing. 

In the sequel, we shall assume that every self-contracted curve 7 : / 1— > M."' is (defined and) 
continuous at the endpoints of /. 

Remark 2.4 (Basic properties). 

(i) Inequality (jl.3p shows that the image of a segment (a, b) by a self-contracted curve 7 lies in 
a ball of radius p := dist(7(a), 7(6)). 
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(ii) A self-contracted curve might not be (left/right) continuous. A simple example is provided 
by the following planar self-contracted curve: 

( if t G (-00,0) 

7(0 = <^ (0,0) ift = 

[ (t,-l) if t £ (0,+oo) 

(iii) If t G (a, b) i-^ 7(t) is a self-contracted curve, then the curve t G (a, b) i-^ 7(0 + b — t) is not 
necessarily self-contracted. 

(iv) Corollary 12.31 reveals that the trajectories of a general gradient system 

at) = -V/(7(t)), 7(0) = xo G 

might fail to be self-contracted curves. Indeed in |16^ page 14] an example of a C°° function 
/ : — > M is given, for which all trajectories of its gradient system are bounded but fail to 
converge. 

(v) In Section [6l we show that whenever / is (quasi) convex, the gradient trajectories are self- 
contracted curves. Thus, bounded self-contracted curves might fail to converge for the strong 
topology in a Hilbert space (see Baillon's example in [2]). 

From now on, we restrict ourselves to the two-dimensional case, and study the asymptotic 
behaviour of self-contracted planar curves. 

3. Horizontal and Vertical directions 

In this section, we introduce a binary- type division of planar segments into horizontal and 
vertical ones. We shall apply this decomposition for segments issued from polygonal line ap- 
proximations of a bounded self-contracted curve. In this section, we derive an upper bound on 
the total length of the vertical segments, while in the next section we shall do the same for the 
total length of the horizontal ones. Combining both results we shall thus obtain an upper bound 
estimation on the total length of a bounded self-contracted curve, establishing Theorem 11.31 

Fix < r < i? and let U{r,R) be the annulus defined in (II. 4p . Let a be a segment 
of U{r, R), not reduced to a point. Denote by p and q the endpoints of a and by m its midpoint. 
Switching p and q is necessary, we can assume that q is closer to the origin O than p, that is 
dist(0,g) < dist(0,p). Let Omq := 6 be the angle between the vectors mO and mq, cf. Fig. [TJ 
Note that 9 € [— f j f ] (by convention, inverse-clockwise angles are positive). 



Lemma 3.1 (Segment length estimate). Let a be a segment ofU{r,R) with endpoints p and q 
such that 9 Then 

2 

length(o-) < |dist(0,p) -dist(0,g)|. 

cos 9 

Proof. Let p be the orthogonal projection of p to the line Om. Using elementary trigonometry 
in the right-angled triangle ppm, we derive 

dist(m,p) = - cos 9 ■ length((T). 

Hence, 

dist(0,p) — dist(0, m) > - cos 9 ■ length(cT). 
Since dist(0,g) < dist(0,m), the conclusion follows. □ 
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Figure 1. 

Fix a € (0, let A G (0, 1) be such that sin a < A < 1, and set r := \R. Denote by 

A := U{XR,R) 

the corresponding annulus of p.4|) . with width equal to AR = (1 — X)R. We now introduce a 
crucial definition in the proof of our main result. 

Definition 3.2 (Classification of the segments). Let a € (0, A € (0, 1) and A as above. A 
nontrivial segment o" of vl is said to be 

• vertical, if 9 lies in (— ^ + a, ^ — q) ; 

• horizontal, pointing in the positive direction, if 9 lies in [— ^, — ^ + a] ; 

• horizontal, pointing in the negative direction, if 9 lies in — a, 

For instance, the segment [p, q] in Fig. [1] points in the negative direction. 

Definition 3.3 (Polygonal approximation). Let 7 : / ^ be a continuous self-contracted 
planar curve converging to the origin O. We consider polygonal approximations {o'i}^^^ of 7 in 
the annulus A, as follows: Let ti < ^2 < • • • < ifc+i be a sequence of points of / with 7(ti) / 
7(ti+i) such that the restriction of 7 to [ti, tfc+i] lies in A. Refining the subdivision if necessary, 
we can further assume that for every i G {1, . . . ,k} the segment at with endpoints pi = 7(ti) 
and Qi = pi+i = 7(tj+i) lies in A and that the length of cij is within any desired precision r/ > 
of the length of 7j[t-,i._^^] (this precision > will be defined at the beginning of Section d] 
and will only depend on a, A and R). Since the function t dist(0,7(t)) is nonincreasing 
(c/. Corollarv 12. 3p . we can further assume that, if cjj := [^4,91] is vertical, then qt is the closest 
point of CTj to the origin. We denote by rrii := ^{pi + qi) the midpoint of cjj. 

Remark 3.4. It is worth noticing that the polygonal approximation of a self-contracted curve 
introduced above is no more a self-contracted curve in general. Nevertheless, one still has 

dist{pi,pi) > dist{pj,pi) when l<i<j<l<k + l, 

which is the property we will use. 

The total length of the vertical segments satisfies the following upper bound. 
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Lemma 3.5 (Total vertical length upper bound). If is a polygonal approximation of j 

in the annulus A, cf. Definition \3.3\ then 

2 

V length(o-i) < AR 

^-^ cos a 

iev 

where the sum is taken over all indices i G V C + corresponding to the vertical 

segments. 



Proof. Let 9i denote the angle between miO and rmqi. Since i € V, it follows that \6i\ < a, 
whence (cos^j)"^ < (cosa)~^. From Lemma |3. II (segment length estimation), we obtain 

2 

y^length((Ji) < ^dist(0,pi) -dist(0,gi). 

COS Ol 

iev iev 
Since dist(0,pj) > dist(0, Qj), for alH G {1, . . . , k] we deduce 

k 

^dist(0,K) -dist(0,g,) < ^dist(0,pi) -dist(0,gi). 

iev i=i 

Now, since qi = Pi+i, the right-hand term is equal to dist(0,pi) — dist(0,gfc), which is less or 
equal to the width Ai? of A. The proof is complete. □ 

4. Length estimate for horizontal directions 

Li this section, we keep the notations and the definitions from the previous section. In 
particular, 

aG(0, |), sina<A<l, A:=U{XR,R) (4.1) 

and {ai}^^i is a polygonal approximation of 7 in the annulus A, cf. Definition 13.31 

We establish an upper bound on the total length of the horizontal segments issued from the 
polygonal approximation of 7. 

Let X & A. The distance from the origin to the half-line passing through x and making 
an angle a > with xO is equal to sin a • dist(0,a;). Thus, the half-line intersects the 
circle S{0, XR) of radius XR centered at the origin at two points. These two points are noted 
7r(x) and vr'(x), with 7r(x) closer to x than tt'{x). Furthermore, the half-line Lx intersects A 
along two segments A^^ and A^, where the endpoints of A^^. agree with x and it{x), and one 
of the endpoints of A^. agrees with 7r'(x). Note that mina;.g^ || 7r(x) — 7r'(x) ||> 6 > 0. The 
half-line Lx extends to a line which bounds a (closed) half-plane Hx containing the origin of the 
plane. Denote by the (closed) half-plane with the same boundary as Hx, not containing the 
origin (see Fig. Ofor an illustration of these notations). The mappings x 1— >■ ir{x) and x 1— > 7r'(x) 
from A to 5(0, XR) are clearly continuous, thus also uniformly continuous. Therefore, there 
exists rj > such that for every pair of points x,y & A which are ry-close from each other (i.e. 
dist(3;,y) < r]), we have 

dist(x, 7r(y)) < dist(x, 7r'(y)), dist(7r(x), 7r(y)) < dist(7r(x), 7r'(x)), 

and dist(7r(x),7r(y)) < dist(7r(x),7r'(y)). (4.2) 



We shall further need the following technical lemmas. 
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Figure 2. 



Lemma 4.1 (Essential disjointness of A^,, Aj,). Let x and y be two distinct points of A such 
that dist{0,x) > dist(0,?/). If y € H^, then the segment Ay lies in Hx- Furthermore, A.y does 
not intersect A^., except possibly at y. 

Proof. Suppose first that y lies in A^^. One easily checks that the angle Ozti{x) increases when 
the point z moves from x to ti{x) along A^;. In particular, the angle OyTT{x) is greater than a. 
Therefore, the segment Ay lies in Hx and meets A^; only at y. 

Suppose now that y ^ A^ and (towards a contradiction) that Ay intersects A^. at z ^ y. Let 
y' denote the intersection point of A^^ with the circle of radius \0y\ centered at the origin. Then, 
the image of Ay by the rotation around the origin taking y to y' does not lie in Hx (the image 
of z should lie in H^). On the other hand, this image agrees with Ayi since the rotation sends 
the ray Oy to Oy'. From the previous discussion, we conclude that the image of Ay is contained 
in Hx- Hence a contradiction. 

Finally, suppose that y ^ A^., Ay n Ax = and Ay intersects A^ at z. li z = Tr'(x), then 

obviously A^^ lies in Hx- Suppose now that z ^ 7r'(x) (thus z ^ vr(y)). Since the angle Oz-K{y) is 

positive while the angle Oztt'{x) is negative, we obtain that 7r'(x)z7r(y) is positive which is not 
possible. The proof is complete. □ 

Lemma 4.2 (Injectivity of vr). Let a := \p,q\ be an horizontal segment of A, with midpoint m, 
pointing in the positive direction. Assume dist(0,g) < dist(0,p). Then, 

(1) the restriction of tt to a is injective; 

(2) if the length of a is at most rj, then the circular arc iT{[p,m]) lies in H^. 

Proof. Let x,y € \p,q] with dist{p,x) < dist(p, y). Since the horizontal segment a points in 
the positive direction, the angle xOy is positive and y lies in Hx. From Lemma l4.ll (Essential 
disjointness of Aa,, Ay), 7r(x) and 7r(y) are distinct (the case y = tt{x) is impossible since it 
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would yield that the angle 7r{x)Ox = yOx is positive, a contradiction). Hence the first part of 
the lemma follows. 

Let X € {m,p]. From above, the midpoint m of o" lies in and the segments A^^ and 
do not intersect each other from Lemma I4.1[ By definition of r/, in view of (j4.2p the points x 
and 7r{x) are closer to 7r(m) than to 7r'(m). Thus, the segment A^,., which does not intersect A^, 
does not intersect A'^ either. That is, A^^ lies either in or in H^. Since the horizontal 
segment a is pointing towards the positive direction, the point x belongs to H^. Therefore, the 
same holds true for the other endpoint tt(x) of A^;. It follows that the circular arc tt{[p, m]) with 
endpoints 7r{p) and 7r(m) is contained in H^. □ 

Lemma 4.3 (Length estimate for horizontal segments). Let a := \p,q] be an horizontal segment 
of A with midpoint m, pointing in the positive direction. Assume dist(0,g) < dist(0,p). Then, 

2 

length((T) < — length(7r([p, m])). 
A 

Proof. The line passing through p and the origin O together with the circle of radius \0p\ 
centered at the origin define a decomposition of the circle of radius \pm\ centered at p into four 
arcs. One of these arcs, denoted by C, contains the point m. Let m' be the endpoint of C lying 
in the circle of radius \Op\ centered at the origin, cf. Fig. [3] below. By construction, 

length((T) = 2 \pm\ = 2 \pm'\. (4-3) 

Since m' is at the same distance from the origin as p, there exists a rotation p centered at the 
origin which takes p to m! . This rotation sends the ray [0,p\ to [0,m'] and preserves distances 
and angles. Therefore, it also sends Aj, to A^/. In particular, the rotation p maps Tr{p) to 7r(m'). 
From Thales' formula, we derive 

|7r(j))7r(m')| |07r(p)| 
\pm'\ \0p\ 

Hence, 

|7r(p)7r(m')| > \ \pm'\. (4.4) 
Since the endpoints of the segment [7r(p), 7r(m')] lie in the arc 7r([pm']), we have 

|7r(p)7r(m')| < length(7r([p, m'])). (4.5) 

When a point x, starting at ml , moves along C, the angle Oxp increases from less than ^ 

to TT. Thus, there exists a unique point m" of C where the angle Om"p is equal to ^ + a. 

By definition of an horizontal segment pointing in the positive direction, the angle Omp lies 
between ^ and ^ + a. Therefore, the point m lies in C between m' and m" , cf. Fig. [3j 

The angles Om"p and Om"Tr{m") are equal to ^ + a and a. Therefore, the ray pm" makes a 
right angle at m" with the line D" passing through m" and Ti{m"). Thus, the line D" is tangent 
to C at m" . This implies that the points O, m and m" lie in the same half-plane delimited by 
the line D' passing through m' and parallel to D" . Since the ray Om' makes an angle less than a 
with D' at m' , the angle between D' and L^' (the half-line passing through m' and making an 

angle a > with m'O) is positive, cf. Fig. [3l Therefore, the points O, m and m" lie in Hm'- 

By applying Lemma [4. II (Essential disjointness of A^,, Ay), successively for x = m' and y = m, 
and for x = p and y = m', we obtain that 7r([p, m']) is contained in tt{[p, m]) from the injectivity 
of the restriction of tt to the segments [p, m'] and [p, m], cf. Lemma l42l Hence, 

length(7r([p, m'])) < length(7r([p, m])). (4.6) 
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Figure 3. 



Putting together the inequahties (|4.3p . (|4.4p . (j4.5p and (j4.6p . we obtain the desired bound. □ 
Let us now consider a polygonal decomposition in A 

of a bounded self-contracted curve 7 converging to O, cf. Definition 13.31 

Lemma 4.4 (Disjointness of TT{pi,mi) and ■7T{pj,mj)). Let di and aj be two distinct horizontal 
segments of a polygonal approximation ofj in A, cf. Definition \3.3l both pointing in the positive 
direction. Then, the images by it of[pi,mi] and [pj,mj] are disjoint. 

Proof. Switching the indices i and j if necessary, we can assume that i < j. 

From Lemma 14.21 the arc 7r([pj,mj]) is contained in H^,. To prove the desired result, it is 
enough to show that 7r([pj,mj]) lies in the complement of {i.e. the interior of HmJ- 

From the definition of a self-contracted curve, the points pj and qj are closer to qi than to pi 
(see Remark 13. 4p . Thus, pj and qj, and so their barycenter ruj, lie in the half-plane delimited 
by the perpendicular bisector of crj. (Notice that this half-plane also contains the origin O, in 
view of Corollarv 12.31 1 The half-line of this bisector with endpoint which makes an acute 
angle with the ray rriiO is noted Drm- Since the horizontal segment crj points in the positive 
direction, its half-bisector makes an angle less or equal to a with ruiO. Thus, lies in 
the half-plane delimited by the perpendicular bisector of crj not containing the origin. 

Now, since the function t 1-^ dist(0,7(t)) is nonincreasing, the points pj and qj, and so their 
barycenter nij, belong to the disk of radius \Oqi\ < \Omi\ centered at the origin. Therefore, the 
points Pj and ruj lie in Hm.., and dist(0,mj) > max{dist(0,pj), dist(0, mj)}. From Lemma HT] 
(Essential disjointness of A^^, A^^), the segments Ap^ and A^^ lie in and do not intersect 

its boundary. Therefore, their endpoints n^pj) and 7r{mj) also lie in the interior Hrm- D 

The total length of the horizontal segments satisfies the following upper bound. 
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Figure 4. 



Lemma 4.5 (Total horizontal length upper bound). If {(yi\\^i is a polygonal approximation of 
7 in the annulus A, cf. Definition \3.3[ then 

J^length(ai) <r^-—AR 

ien 

where Ti is the set of indices corresponding to the horizontal segments. 

Proof. From Lemma 14.31 (Length esthnate for horizontal segments), the sum of the lengths of 
the horizontal segments pointing in the positive direction satisfies 

^ length(cri) < ^ length (7r([pj,mj])), 

where 7^+ is the set of indices corresponding to the horizontal segments pointing in the positive 
direction. Since the arcs 7r([pj,mj]) of the circle 5(0, Ai?) are disjoint, cf. Lemma 14.41 we have 

length(7r([pi,mi])) < 2t:XR. 

Analogous arguments lead to a similar bound for the sum of the lengths of the horizontal 
segments ai pointing in the negative direction. Recalling that AR = (1 — X)R the result follows. 

□ 

5. Proof of the main result 
In order to prove our main theorem (cf. Theorem 1 1.3p . we shall first need the following result. 

Proposition 5.1 (Length estimate in the annulus A). Every continuous self- contracted planar 
curve 7 converging to the origin O satisfies 

length(7 H A) < (Svr + 2) AR. 

Proof. Consider a decomposition of 7 into segments cjj as in Definition 13.31 (refining a subdi- 
vision does not decrease the sum). From Lemma 13.51 (Total vertical length upper bound) and 
Lemma 113] (Total horizontal length upper bound), the length L of the polygonal line C7ia"2 • • • cr^ 



12 



A. DANIILIDIS, O. LEY, AND S. SABOURAU 



satisfies 

L = ^ length ((Ti ) + ^ lengtli (di ) 
\ 1 — A cos a J 

By taking the supremum of L over all such decompositions respecting the annulus A, we derive 
the same upper bound for the length of 7 fl ^. Finally, by letting a and A go to zero, we obtain 

length(7 r\ A) < (Svr + 2) Ai?. 

□ 

Now, we can derive our main result. 

Proof of Theorem \1.3l From Corollary 12.31 the bounded continuous self-contracted curve 7 con- 
verges to a point. Using a translation if necessary, we can assume that this point agrees with 
the origin O of . 

Let t- = inf /. Denote by 7(t-) the limit of ^{t) when t goes to t_, (c/. Proposition 12.21 
(Existence of left/right limits)). Set Rq = dist(0, 7(t_)). For i G N, let Ai be the planar 
annulus centered at the origin with outer radius Ri and inner radius where Ri+i = XRi, 

with A G (0, 1) given in (j4.ip . From Proposition 15.11 (Length estimate in the annulus), we have 

length(7 n Ai) < {Sir + 2) ARi (5.1) 

where ARi is the width of Ai. Since A < 1, the sequence Ri goes to zero and the sum of the width 
of the disjoint annulus Ai is equal to Rq. Thus, taking the sum of the above inequalities (j5.ip 
for i G N we obtain the desired bound 

length(7) < (Svr + 2) dist(0, 7(t-)). 

The proof is complete. □ 

6. Gradient and subgradient systems, and convex foliations 

In this section, we apply Theorem 11.31 to derive length estimates, first for orbits of dynamical 
systems of gradient or subgradient type, then for orbits orthogonal to a convex foliation. The key 
fact is to observe that in some interesting particular cases (for instance, / convex or quasiconvex) 
these curves are self-contracted. Recall however that this is not the case for gradient dynamical 
systems defined by a general C°° function, as already observed in Remark 12.41 (iv). 

6.1. Gradient dynamical system — quasiconvex case. Let / : M" ^ M be a C'^ function 
(/c > 1), xo G and consider the Gradient Dynamical System 

= -V/(7(t)), t>0 
7(0) = xo G M". ^ ' ' 

It follows directly from the standard theory of Ordinary Differential Equations (see [12], for 
example) that the system ()6.ip admits a solution (trajectory) 7 : / 1 — > M", where I C [0, -|-oo), 
which is a curve of class C^~^. Note that the case k = 1 corresponds to mere continuity of 7, 
while for k > 1 (or more generally, if / is assumed C^'^, that is, V/ is Lipschitz continuous), 
the trajectory 7 is unique. In the sequel, we shall always consider maximal solutions of (j6.ip . 
that is, for which / = [0, T), where T > is the maximal time such that 7 is defined in [0, T). 
We shall refer to them as orbits of the gradient flow of /. 
We will also need the following definition. 
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Definition 6.1 (Convex, quasiconvex and coercive functions). A function / : M" — > M is called 
convex (respectively, quasiconvex) if its epigraph 

epi/:={(x,y) eM"+^ | f{x)<y} 

is a convex subset of x R (respectively, if for every y € M the sublevel set {x G M" | f{x) < y} 
is a convex subset of R"). A function / is called coercive (or proper), if its level sets are bounded, 
or equivalently, if 

lim f{x) = +00. (6.2) 

>+oo 

It is straightforward to see that whenever / is coercive the corresponding flow orbits are 
bounded curves and therefore / = [0, +00) (the trajectories are defined for all t > 0). Notice in 
particular that the function 

t^filit)), t€[0,+oo) 
is a natural Lyapunov function for the orbits of the flow, i.e. it is nonincreasing along the 
trajectories. Moreover, unless 7 meets a critical point {i.e. V/(7(t*)) = for some € [0, +00)), 
the function defined in (j6.ip is decreasing and 7 is injective. 

Let us finally recall {e.g. [9, Theorem 2.1]) that a (differentiable) function / : R" — > R is 
quasiconvex if and only if for every x, y E R" the following holds: 

{Vf{x),y-x)>0^ f{y)>f{x). (6.3) 

We are now ready to establish the following result. 

Proposition 6.2 (Quasiconvex orbits are self-contracted curves). The orbits of the gradient 
flow of a quasiconvex C^'^ function are self- contracted curves. 

Proof. Let 7 : / 1 — > R" be an orbit of the gradient flow of /. Let < t < ti be in / and consider 
the function 

9{t) = l\h{t) - 7{t,)\\\ tel. 

In view of ()6.ip . we easily deduce that 

9'{t) = {Vf{^{t)),^{t^)-^{t)). 

If g'{t) > for some t € [0, ti), then the quasiconvexity of / would imply that f{'~f{ti)) > f{'y{t)) 
(see dSiSI)), which in view of i^Ji) would yield that -f{t') = 7(ti) for all t' £ [t, ti] and Vf{^{t)) = 
0, a contradiction. Thus, g is nonincreasing in the interval [0, ti], which proves the assertion. □ 

The following corollary is a straightforward consequence of the previous proposition and 
Theorem 11.31 (Main result). 

Corollary 6.3 (Orbits of a gradient quasiconvex flow). Let f : M? ^ W be a coercive C^'^ 
quasiconvex function. Then, for every xq € M^, the orbit of the gradient flow (|6.ip converges 
and has finite length. 

6.2. Subgradient dynamical systems - convex case. A convex function is a particular 
case of a quasiconvex function. Therefore, Corollary 16.31 implies that the orbits of the gradient 
flow of C^'^ convex functions are of finite length. It is well-known ([6J) that in the case of a 
(nonsmooth) convex function / : M" — > M (or more generally, for a semiconvex function [lOj), the 
gradient system (j6.ip can be generalized to the following differential inclusion, called Subgradient 
Dynamical System 

r 7(0 e -df{l{t)) a.e. t € [0, +00), 

\ 7(0) = xo G R", ^ ^ 
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where df is the set of the subgradients (subdifFerential) of /. If / : M" ^ M is convex, then this 
latter set is defined as 

df{x) = {peW\ f{y) > f{x) + {p,y-x),yye W} for all xeW. 

The above formula defines always a nonempty convex compact subset of M", and reduces to 
{V/(x)} whenever / is differentiable at x, c/. [8]. It is also known that (j6.4p has a unique 
absolutely continuous solution 7 : [0, +00) — M", that is, the derivative ^{t) = ^7(t) exists 
almost everywhere and for every < ti <t2, 

7(t2) = 7(ti)+ r i{t)dt. 
Jti 

The analogous of Proposition 16.21 holds true. 

Proposition 6.4 (Subgradient convex flow orbits are self-contracted curves). Let f : — M 

be a convex continuous function. Then, for every xq £ M", the trajectory 7 of the subgradient 
system (j6.4p is a self-contracted curve. 

Proof. We give a sketch of proof for the reader convenience (we refer to [6j for details). It is 
easy to prove that t G [0, +00) f{'y{t)) is convex and that for almost all i > we have 

|/(7W) = -||7(i)lP<0. 

Therefore t filit)) is nonincreasing and j{t) G {/ < /(xq)} is bounded. Moreover, for all 
ti > and for almost all t G [0, ti] 

\ 1"^^*) " = (7W,7W - l{ti)) < /(7(ti)) - f{l{t)) < 0. 

This implies that t G [0, ti) 1— > Il7(i) — 7(ii)|P is nonincreasing yielding that 7 is a self-contracted 
curve. □ 

When n = 2, we have the following generalization of Theorem ll.il 

Corollary 6.5 (Orbits of a subgradient convex flow). Let f ■.R'^ be a convex continuous 

function which admits a minimum. Then, for every xq G M^, the orbit of the gradient flow (j6.4p 
converges and has finite length. 

6.3. Trajectories orthogonal to a convex foliation. In this section we consider orbits that 
are "orthogonal" to a convex foliation. Let us introduce the definition of the latter. (For any 
subset C C M^, int C denotes the interior of C and dC its boundary.) 
Let {C'alogjo^A] (where A > 0) be a family of subsets of such that 

(i) For all a G [0,^], Ca is convex compact 

(ii) If a > a', then Ca C int Ca', (6-5) 

(iii) For every x G Cq \ int Ca, there exists a unique a G [0, A] such that x G dCa- 

We shall refer to the above as a foliation made up of convex surfaces. We shall now define a notion 
of "orthogonality" for an orbit 7 with respect to this foliation. To this end, let T G (0, +00] and 
7 : [0, T) — > be an absolutely continuous curve. We say that the curve 7 is "orthogonal" to 
the foliation defined in ()6.5p if the following conditions hold: 

(i) for every t G [0, T), there exists a G [0,A] such that j{t) G dCa, 

(u) for almost ah t G (0, T), if j{t) G dCa, then for ah x G Ca, {i{t),x - 7(t)) > 0, (6.6) 
(iii) ift'>t and j{t) G Ca, then j{t') G Ca- 
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Condition (ii) in (j6.6p is a nonsmooth generalization of orthogonality: if dCa is smooth at ^{t) 
and 7 is differentiable at t then ^{t) is orthogonal to the tangent space of dCa at 7(i). Further, 
condition (iii) guarantees that the curve 7(t) enters into each convex set of the foliation. In this 
context, one has the following result. 

Proposition 6.6 (Orbits orthogonal to a convex foliation). The curve ^ is a hounded self- 
contracted curve, thus, of hounded length. 

Proof. The curve 7 is clearly bounded. Let < ti < T. Then, for almost all t € [0, ti], we have 

\j^Mt)-i{ti)\? = m),i{t)-i{ti)). (6.7) 

By (16.61) (i), we have 7(t) £ dCa for some a. By (16. 6p (iii) and since ti > t, we also have 
7(^1) S Ca- Therefore, (16. 6p (ii) implies that the right-hand side of ()6.7p is nonpositive. It 
follows that t G [0,ti] i—f \\^{t) — 7(^1) |P is nonincreasing and 7 is self-contracted. Applying 
Theorem 11.31 finishes the proof. □ 

Remark 6.7. (i) The sublevel sets of a continuous quasiconvex function need not define a 
convex foliation. Indeed, consider the quasiconvex function / : [—2, 2] — M given by 

{X, if - 2 < X < 0, 
0, if < X < 1, 
X - 1, if 1 < X < 2. 

Then the sublevel sets of / do not define a foliation on [—2, 2] C M since property (iii) of (j6.5p 
fails at the level set [/ = 0] . This drawback appears whenever the level sets of such functions have 
"flat" parts outside the set of their global minimizers. Actually, it follows from ^ Theorem 3.1] 
that the sublevel sets of a continuous quasiconvex coercive function / define a convex foliation 
if and only if the function is semi- strictly quasiconvex. (We refer to [9] for the exact definition 
and basic properties of semi-strictly quasiconvex functions.) 

(ii) Let / : — >■ M be a coercive C^'^ quasiconvex function and 7 : [0, +00) — > be the 
solution of (16. ip . Let Xqo be the limit of ^{t) as t ^ +00 and assume that / has no critical 
point in {/(xqo) < / < /(^^o)}- Then, it is not difficult to see that {/ < o}aG[/(xoo),/(a;o)] ^ 
family of convex compact subsets which satisfies (|6.5p (in fact, / is semi-strictly quasiconvex 
in [/(xoo), /(xo)]) and 7 satisfies (j6.6p . 

(iii) Despite the first remark. Proposition 16.61 can be used to obtain the result of Corollary 16.31 
without the extra assumption made in the second remark. The reason is that the trajectory 
of the gradient flow will not pass through the flat parts of / anyway (if it does, then it stops 
there). We leave the technical details to the reader. 



7. Two COUNTER-EXAMPLES 

7.1. Absence of Convexity. The second conclusion of Corollary 16.31 fails if / is not quasicon- 
vex, even when the function is C°° and has a unique critical point at its global minimum. Let 
us give an explicit example below: 

Define a function / : — > M in polar coordinates as 

f^r,e)=e-'/^{l+r + sm{}.+e)) 
with f{0) = 0. The graph of / in the plane 6 = looks like the graph of Fig. [5l 
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Figure 5. 



One can check that / is smooth, positive away from O, with no critical point except at the 
origin. The gradient trajectory of / issued from the point {r,9) = ((^)~^,0) remains close to 
the spiral given by 

r = (f 

e = -t 

where t runs over [0, oo). Therefore, it converges to the origin and its length is infinite. 

7.2. Thom conjecture fails for convex functions. Let / : ^ M be a convex continu- 
ous function which admits a minimum. Then, Corollary 16.51 guarantees that the orbits of the 
gradient flow of / have finite length (thus, a fortiori, are converging to the global minimum of 
/). However, it may happen that each orbit turns around its limit infinitely many times. In 
particular the corresponding statement of the Thom conjecture fails in the convex case. 

We construct below a counter-example using a technique due to D. Torralba [18] which allows 
us to build a convex function with prescribed level-sets given by a sequence of nested convex 
sets. Let us recall his result. 

For any convex set C C M", the support function of C is defined as 6cix*) = sup2.g(;7(x, x*) 
for all X* E M". Let {Ck}keN be a decreasing sequence of convex compact subsets of such 
that Ck+i C intCfc. Set 

Kk = max - — — — -. 

11^*11=1 ^Ck{x*) - 5c^+^{x*) 

Then Torralba's theorem [18] asserts that for every real sequence {AfcjfcgN satisfying 

< Kk{\k - Afc+i) < Afc_i - Afe for every k>l, (7.1) 

there exists a continuous convex function / such that for every A; G N, {/ < A^} = Ck- Moreover, 
Afc converges to min/ and, for any A; > and A G [Afc+i, Afc], we have 

A — A/c+i \ ^ ( Afc — A 



{/ < A} = ft^ \Ck+{ ■ ' ■ Ck+i (7.2) 

(i.e., the level-sets of / are convex interpolations of the two nearest prescribed level-sets). 

Ste'p 1. Constructing a first piece of trajectory. Consider the finite decreasing sequence of convex 
sets Cq = B{0,1), Ci = B{0, 0.9), Cs = C3 = 5(0,0.6) and d = B{0, 1/2) where E is a 
compact set bounded by an ellipse (see Fig. [6]). It is easy to find a sequence {Xk} which satisfies 
(|7.ip : set K = mayi{Ki, K2, K3, K4} + 1 > 1 (since C C int C implies 6c > ^c), take Aq = 1, 
Ai = 1/2 and 



Afc — Afc+i — -^(Ao — Ai). 



(7.3) 




Figure 6. 



We then obtain a positive function /o : Co — > M with argmin /o := {/o = niin /o} = C4. Consider 
the subgradient trajectory 70 starting from the point Aq of Co (see Fig. [6]). It reaches S 9C4. 
From ()7.2p this trajectory is radial (pointing towards the origin) between dCo = {/o = Aq} and 
dCi = {/o = Ai} and between dC^ = {/o = A3} and dC^ = {/o = A4}. Due to the presence 
of the ehipse C2, the trajectory deflects with an angle 9 := A0OA4 > in the clockwise direction. 

Step 2. Construction of the function from the previous step. Consider the transformation 
T = r o h, where h is the homothety of center O and coefficient 1/2 and r is the rotation of 
center O and angle 9. We define, for all A; G N and k £ {0, 1, 2, 3} 

Cfc = T'^'/'^I (C^) where [k/4] is the integer part of k/A and k = k (modulo 4) 
(see Fig. [7] for the ffist steps of the construction). 

The sequence of convex sets {C^} satisfies the assumptions of Torralba's theorem and we can 
define a sequence {A^} as in ()7.3p which satisfies (jT.ip (note that {Kk} is 4-periodic since, for 
all convex set C C and x* G M?, Sq-(^c)i^*) = ^^cir'^ix*))). We obtain a convex continuous 
function / : Co — > M^" with argmin / = {O}. The trajectory starting from the top of Cq spirals 
around the origin and converges to O (see Fig. [7]where the beginning of the trajectory is drawn 
with a deflection of 30). 

Step 3. Smoothing f. Actually, the function / built above is C°° except at the origin and at 
the boundaries dCk- It is possible to smooth out / in order to obtain a function which is C°° 
everywhere except at the origin and C"^ at the origin (for any fixed m > 1). The smoothing 
procedure is quite involved from a technical point of view and is omitted. We refer the interested 
reader to [H Section 4.3] where such a smoothing is realized (in a different context). This 
procedure does not modify significantly neither the function nor its gradient trajectories (i.e. 
they remain a spiral). This concludes the construction. 
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Figure 7. 
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